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The  purpose  of  this  thesis  was  to  develop  a  computer 
program  to  model  the  motion  and  radar  beam  characteristics 
of  different  configurations  for  the  distributed  sparse 
array.  The  effect  of  the  two-body  motion  on  the  beam 
quality  of  the  array  was  analyzed.  Two  groups  of  arrays, 
planar  and  three-dimensional  were  considered.  The  planar 
phased  arrays  were  rectangles  and  disks,  and  the  other 
arrays  included  cones  and  spheres.  The  number  of  emitters 
in  the  configurations  ranged  from  five  to  two  hundred. 

The  beam  quality  parameter  was  the  half-power 
beamwldth  as  determined  in  discrete  directions  throughout 
the  hemisphere  below  the  array.  An  array  was  considered 
feasible  If  the  beamwidths  degraded  only  slightly  during  a 
orbit. 

The  most  useful  configuration  was  found  to  be  a  spher 
of  randomly  spaced  emitters  because  it  provided  narrow 
beamwidths  in  all  directions  below  the  array.  Because  of 
the  motion  of  the  array,  the  beamwidths  changed  in  a 
predictable,  periodic  manner.  The  half-power  beamwidths 
actually  improved  in  certain  directions  and  only  slightly 
degraded  in  others.  The  sphere  also  required  the  least 
number  of  emitters  and  could  be  used  at  any  altitude. 


EFFECT  OF  TWO-BODY  MOTION  ON  RADAR  BEAM  QUALITY  FOR  VARIOUS 


DISTRIBUTED  SPARSE  ARRAY  CONFIGURATIONS 


L-l.  Introduction 


Background 

In  1586  the  USAF  unveiled  a  study,  Project  Forecast  II, 
that  established  chosen  technologies  and  systems  concepts  as 
initiatives  for  the  Air  Force  Systems  Command  to  pursue  and 
for  operational  commands  to  support.  The  study  identified 
thirty-nine  technologies  and  thirty-one  advanced  systems 
concepts  which  would  "revolutionize  the  way  the  Air  Force 
carries  out  its  mission  in  the  twenty-first  century, 

A  third  of  the  systems  concepts  listed  were  related  to 
spacecraft  or  space  missions.  This  thesis  is  concerned  with 
the  concept  of  a  "distributed  sparse  array  of  spacecraft." 
The  idea  "involves  placing  large  phased  arrays  in  space  with 
major  components  of  the  array  not  rigidly  connected  to  each 
other"  (2:47,49  )  . 

Phased  array  radars  are  ideal  for  search  and  tracking 
of  targets  over  large  areas  because  of  the  ability  to 
electronically  steer  the  radar  beam.  Hundreds  of  targets 
can  be  tracked  almost  simultaneously  by  moving  the  beam  from 
target  to  target  in  a  matter  of  microseconds.  An  example  of 
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such  capabilities  is  the  PAVE  PAWS  _adar  system.  A  single 

array  can  span  120  degrees  in  azimuth  and  follow  a  large 
number  of  objects  almost  simultaneously  by  electronically 

shifting  its  beam  from  one  target  to  another  within  a  few 
millionths  of  a  second  (1:95). 

The  distributed  sparse  array  radar  network  could  be 
used  for  detection  and  early  warning  of  bomber  or  cruise 
missile  attacks.  Cruise  missiles  are  a  threat  because  they 
can  fly  close  to  the  surface  under  radar  beams  and  around 
defenses.  Look-down  sensors,  particularly  space-based 
radars  provide  feasible  means  of  detecting  the  missiles. 

The  radars  could  expand  the  earth  coverage  for  better 
warning  and  tracking  (3:78). 

The  distributed  sparse  array  also  introduces  a  new 
degree  of  survivability.  The  network  would  consist  of 
several  nonsophist icated  and  relatively  inexpensive 
satellites  which  would  reduce  the  reliance  on  a  few 
extraord  inar i ly  capable,  expensive  and  vulnerable 
satellites.  "It  therefore  may  be  possible  to  create  a 
phased  array  device  that  we  can  place  into  space  and  enhance 
simply  by  adding  more  relatively  Inexpensive  elements 
whenever  the  threat  increases  and  budget  pressures  permit" 
(2:49). 

The  idea  of  a  large  phased  array  in  space  is  appealing 
but  also  very  complicated.  For  a  distributed  pnased  array, 
each  element  would  be  In  its  own  orbit  and  therefore  the 
array  would  tend  to  drift  apart  with  time.  Placing  a  large 
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planar  array  in  space  and  maintaining  the  elements  to  fixed 
relative  positions  would  require  enormous  amounts  of  fuel 
onboard.  A  planar  array  would  therefore  not  be  cost 
effective  if  employed  as  a  distributed  sparse  array. 

The  objective  of  this  thesis  was  to  analyze  possible 
configurations  and  determine  how  the  radar  beam  changed 
during  an  array's  orbit.  The  array  needed  to  form  a  beam 
with  an  adequate  beamwldth  in  any  direction  from  the  array 
at  any  time  in  the  orbit. 

AppppacJb 

To  analyze  the  effect  of  the  relative  motion  of  the 
satellites  on  the  beam  quality,  a  computer  program  was 
developed.  The  program  used  the  locations  of  the  emitters, 
determined  the  phase  differences,  and  calculated  the 
magnitude  of  the  electric  field  in  a  given  direction.  In 
incremental  angular  steps,  the  field  in  various  directions 
was  evaluated  and  used  to  determine  the  half-power  beamwldth 
( HP BW ) .  The  HPBW  was  then  calculated  throughout  the 
hemisphere  below  the  nrray  each  time  the  array  was  moved  in 
its  orbit.  Plots  of  the  HPBW  values  were  analyzed  to  find 
configurations  and  initial  conditions  which  caused  little 
change  in  HPBW  values  during  the  orbit  of  the  array. 
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1 1  ■  Qi-bll  MgUon  Theory 


Background 

The  configuration  of  the  distributed  sparse  phased 
array  is  a  network  of  several  satellites  within  a  region  of 
space.  The  satellites  are  in  separate  orbits  and  therefore 
the  relative  positions  to  the  center  of  the  array  will 
change  as  the  satellites  progress  around  in  their  orbits. 

The  electric  field  of  the  array  depends  on  the  relative 
positions  of  the  emitters.  For  the  array  configuration  to 
remain  effective,  the  electric  field  pattern  needs  to  remain 
relatively  unchanged  throughout  the  orbit. 

To  determine  the  effect  of  the  motion,  the  relative 
locations  of  the  many  satellites  at  any  time  in  their  orbits 
needed  to  be  calculated.  For  the  relative  positions  of  the 
array,  the  motion  of  the  satellites  was  with  respect  to  the 
center  of  the  array.  The  center  of  the  array  configuration 
was  assumed  to  be  in  a  circular  orbit.  The  orbital  analysis 
was  based  on  the  Clohessy-Wi ltsh ire  equations  (10:1). 

Peyelpcment 

The  motion  of  the  array  was  developed  by  first 
considering  its  kinematics.  The  array  elements  were 
considered  point  masses  about  the  central  emitter  which  was 
assumed  to  be  a  point  mass  in  motion  in  a  rotating  reference 
frame  fixed  to  the  earth.  The  vector  positions  of  each 
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satellite  with  respect  to  the  center  of  the  earth  were 
written  as  the  sum  of  vectors  to  the  center  of  the  array  and 
then  to  the  individual  satellites: 

r  =  r  +  p  ( 1 ) 

V  O  V 

The  distances  between  the  points  were  small  relative 

_  A  A  A 

to  so  that  the  displacements  in  the  e^,  e e directions 
were  small.  The  relative  position  vectors  p  were  then 
broken  into  radial,  tangential  and  vertical  components. 
Figure  1  shows  the  r.  ,  ea,  e  coordinate  system  used  for  a 
satellite  in  the  array.  The  center  of  the  coordinate  system 
was  the  reference  satellite,  and  the  vector  to  point  2  was: 

r  =  (r  ♦  6r)e  +  It  66)efl  ♦  6ze  (2) 

2  o  r  o  Q  2 


Figure  1.  Orbital  Coordinate  Frame 
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The  e  ,  e~,  e  axis  frame  rotates  about  the  earth  with 
V  tr  Z 

the  orbital  angular  velocity  co1,  =  coez . 

Taking  two  time  derivatives  of  equation  (2)  would  give 
the  inertial  acceleration  of  each  satellite  about  the 
center.  Equation  (3)  was  the  notation  for  the  absolute  rate 
of  change  of  a  rotating  system  where  (r)^  was  the  rate  of 
change  relative  to  the  rotating  system  (4:48). 

r=(r)r+wxr  (3) 

The  velocity  of  a  satellite  was  expressed  as: 

—  •  •  .  A  •  *■  »  A  «  A  —  — 
r  =  fr  +6rje  +  r  6 0ert  +  t  60e^  +  6ze  +  to  x  r  (4) 

2  v  J  X  o  3  o  &  Z 


where 


to  x  r  =  oo  (ro  +<5r)ee  -  W(ro60)v 


By  taking  another  time  derivative,  the  inertial  acceleration 
was  expressed  as: 

r  =  f<5r  +  tor  <5©  -  cor  66  -  tor  dele 

2  t  o  o  o  '  X 

+  (r<58  ■*•1^60  +  co(r^+6r)  +  tor^  +  <o6r)eQ  +  6zez  +  to  x  r  (5) 
where 


u>  x  r 


toro6e)e0 


For  satellites  in  two-body,  circular  orbits,  the 
angular  and  radial  velocities  were  constant. 

Therefore : 

*  • 
oo  =  0  and  r  =0 
o 

The  mean  orbit  period  was  IP  =  2  n(a*/>j)i/2 
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The  angular  velocity  u>  is  the  rate  the  central 
satellite  completes  an  orbit.  Because  w  =  •  the  angular 

velocity  equals  the  mean  motion  (w  =  n) .  These  relations 
reduced  equation  (5)  to 

r_=[V5r  -  2nr  66  -n*(r  +6r)"]e 

2  •  o  o  J  r 

+  fr  66  +  2n<5r  -  n2r  6©"]e.-  +  6ze  (6) 

L  o  o  J  6  Z 


Assuming  the  array  was  in  a  two-body  orbit  then  equation 
(7)  expresses  the  gravitational  acceleration  (10:2). 


a 
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(7) 


The  magnitude  of  r  =  r  +  p  was: 

2  o 

r2  =  Oo  +  2iO*~P  +  72 


Therefore, 


2r  op  z 

From  equation  (8),  c  was  defined  as  £  =  — — —  +  J-^-| 

r2  t  ro^ 

O 

—  3/2 

The  term  (1+c)  was  expanded  with  the  binomial 

expansion.  Keeping  first  order  terms  of  p  yielded: 


(9) 


Now  r  was  put  in  terms  of  the  e  ,  e_,  e  coordinate  system, 
2  r  “  2 


:  *  a  (r7  +  2r  6r  +  6r2  +  r z6©2  +6Z2") 
2  v  o  O  O  ' 


a -*-•/* 


(10) 
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pi  UU  J  11*1,1 


The  6r,  66,  and  6 z  were  small  with  respect  to  r  so 

O 

2  2  2 

that  6x  ,  66  ,  and  6z  were  approximately  equal  to  zero. 
Equation  (10)  was  rewritten  as: 

r  s  s  r  3(l  +  c)  3</z  where  e  %  (11) 

2  o  v  /  r 

o 

By  using  the  binomial  expansion  as  before,  equation  (11) 
was  reduced  to: 

r's*  r~a(l  -  36r/r  )  (12) 

2  O'-  r  O' 

This  expression  for  r was  used  in  the  denominator  of 
equation  (7)  to  give: 


V6r)*t 


r©6See 


+  <5ze 


(13) 


The  squares  of  the  6  terms  were  3:  0,  and  for  a  circular 
orbit: 

(  1 1/2 

a  =  r  and  n  =  — 1 

°  l  r*J 

O 

Equation  (13)  was  then  simplified  to: 


*  -nZ[(ro"2<5Oer 


zJ6e6 


+  <5ze_ 


(14) 


Assuming  the  array  underwent  only  two-body  motion, 
then  the  gravitational  acceleration  and  the  kinematic 
accelerations  were  equal. 
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r  gave : 


Setting  a  = 

9 

(<5r  -  2nr ^69  -  3nZ<5r)er 

+  £r  JS6  +  2n6r)e0  +  (6z  +  nZ6z)ez  =  0  (15) 


Equation  (15)  shows  the  motion  in  the  e  direction  is 

z 

uncoupled  from  the  e  ,  and  directions.  The  solution  for 
the  6z  displacement  was  found  to  be  a  simple  harmonic 
oscillator.  The  initial  conditions  when  t  =  0  were: 


6z(0)  =  <5z  and  <5z  ( 0 )  =  6z 

o  o 

There  fore, 

<5z 

6z  ( t )  =  6z  cos(nt)  +  — -  sin(nt)  (16) 

o  n 

<5z  ( t )  =  ~n6z  sin(nt)  +  <$z  cos(nt)  (17) 

o  o 

The  eQ  components  were  integrated  with  the  initial 
conditions  of  <5©(0)  =  66  and  66(0)  =  66  to  give: 

O  O 


r  66  { t)  =  2nf  6r  -  drtt)}  +  r  66 

©  V  o  J  o  o 


(18) 


Equation  (18)  was  used  in  the  component  to  give 

<5r  -  4nz6r  +  nz6r  -  2nr  66  =  0  (19) 

o  o  o 


With  the  initial  conditions,  6 r(0)  =  6r  and  6r ( 0 )  =  6 r 

O  O 

the  homogeneous  and  particular  solutions  were: 

6r(t)|h  =  A*cos(nt)  +  B'sin(nt)  (20a) 

6r(t)|  =  46r  +  £/  'Jr  66  (20b) 

1  p  o  v  o  o 
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The  complete  solution  for  the  er  direction  was  then 


6r(t)  =  (-36ro  -  2/^( ro6©o))cos ( nt ) 

6K  2 

+ -  sin(nt)  +  46z  +  /  r  66  (21) 

n  o  r»  o  o 

6r(t)  =  (3n6ro  +  2rQ(5£o)sin(nt)  +  dr^cosint)  (22) 

Using  the  solution  for  6r ( t )  in  equation  (18)  gave  the 

solution  for  r  66( t): 

© 

r  6e(  t)  =  -3r  66  -  6nr  6r  +  |6n6r  +  4r  66  jcos(nt) 

©  OO  00^0  ooj 

-  (26ro) s  i n (  nt  )  (23) 

Integrating  equation  (23)  gave  the  so.ution  for  rQ60(t). 


r  66  ( t )  =  r  66  -  |3r  <56 

o  0  0^0 


66  +  6n6r 

o  o 


It  +  [(>n6z  +  — r  66  lsi 

J  t  o  n  o  oj 


sin(nt ) 


S4 


cos ( nt )  - - 

n 


Equations  16,  17  and  21  -  24  represent  the  equations  of 

motion  in  the  e  ,  e~,  and  e  coordinate  system  for  a 
r  t?  z 

satellite  in  an  orbit  near  the  center  of  the  array.  For  use 
with  several  satellites  the  equations  were  written  in  matrix 
form  in  equation  (25)  (10:2-6). 
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'  <5r  { t ) 
6R  ( t )  =  r  69  ( t ) 

V  O 

6z  ( t ) 


<5r(t) 

and  6V  (t)  =  -r  6©(t) 

l  O 

<5z  ( t ) 


<5R  (t)  =  *  6r(t=0)  +  ft  6v(t=0)  (26a) 

v  rr  rv 

6V  (t)  =  *  6r  ( t  =  0  )  +  ft  <5v  ( t  =  0  )  (26b) 

v.  vr  vv 


The  ft  matrices  were  defined  in  equations  (27a-d)  where 

XX 

'I'  =  n*',  and  t  was  in  seconds  (10:6). 


{  = 

rr 


(  A  -3cos*) 
6  (  s  i  n*  -  *) 
0 


0  * 
0 

cos* 


(  27a  ) 


/  (sin*)  /  fl  -  cos*') 

n  n  v  J 

2/n(cos*  -  l)  1/n(4sin*  -  3*) 


/nSin» 


(27b) 


3n(sin*) 

6n  [cos*  -  l) 


-n  (  s  i  n*) 


(27c) 


fcos* 
-2 ( sin*) 


2s  i  n* 

-3  +  4cos* 


cos* 


(  27d  ) 


The  only  variables  in  these  equations  were  the  initial 

positions  and  velocities,  orbit  altitude  and  amount  of  time 

the  array  moved.  Knowing  these  parameters,  the  matrices 

<5R  (t)  and  6V  ( t )  gave  the  displacements  in  the  e  ,  ea,  and 
v  v  r  B 


e  directions  after  the  satellites  had  moved  t  seconds 
z 
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III..  Array  Antenna  Theory 


B,a.c.Kqi.flund 

Phased  array  radars  are  ideal  for  search  and  tracking 
of  targets  over  large  areas  because  of  the  ability  to 
electronically  steer  the  radar  beam.  Hundreds  of  targets 
can  be  tracked  almost  simultaneously  by  moving  the  beam  from 
target  to  target  in  a  matter  of  microseconds.  A  phased 
array  radar  works  by  a  group  of  identical  emitters  each 
radiating  its  energy  from  the  array.  Depending  on  the  shape 
of  the  antenna,  the  radiation  forms  a  narrow,  pencil-like 
beam,  suited  for  tracking  or  a  fan-like  beam,  best  for 
searching  broad  areas  (1:94). 

When  all  the  signals  leave  the  array  in  phase,  they 
will  add  in  phase  at  any  point  along  a  line  perpendicular  to 
the  plane  of  the  array.  The  signals  constructively 
interfere  (add)  along  the  array's  boresight,  or 
perpendicular  axis,  and  within  a  small  angle  to  each  side. 

At  greater  angles  to  the  boresight,  individual  signals  from 
different  radiating  elements  must  travel  different  distances 
to  reach  a  target.  As  a  result  the  relative  phases  are 
changed  and  the  signals  interfere  destructively.  Thus, 
outside  the  narrow  cone,  centered  on  the  arra/'s  boresight, 
targets  produce  no  detectable  return.  Because  of  the 
character  1st ics  of  interference  patterns,  the  width  of  that 


cone  is  directly  proportional  to  the  operating  wavelength 
and  inversely  proportional  to  the  size  of  the  array  (1:96) 

Beam  steering  is  accomplished  when  the  signals  from 
each  of  the  radiating  elements  are  delayed  electronically  by 
amounts  that  increase  steadily  across  the  array.  Each  delay 
causes  a  signal  to  lag  a  fraction  of  a  wavelength  behind  the 
signal  from  the  adjacent  element.  The  zone  in  which  the 
signals  add  up  in  phase  to  produce  a  return  signal  lies  not 
down  the  boresight  of  the  antenna,  but  off  to  the  side  in 
the  direction  of  Increasing  phase  delay.  The  angle  of  the 
beam  reflects  the  magnitude  of  the  phase  shift,  the  size  of 
the  array  and  the  wavelength  of  the  signals.  The  beam  has 
the  form  of  a  slender  cone  surrounded  by  regions  of 
destructive  Interference  (1:96). 

Antenna  EflUeu 

The  distributed  sparse  array  network  is  a  general  form 
of  a  phased  array  radar.  To  determine  the  beam 
characteristics  of  a  phased  array  radar,  the  radiation 
patterns  from  all  the  emitters  are  combined  to  give  the 
overall  antenna  pattern.  The  antenna  pattern  is  then  used 
to  show  the  direction  and  relative  magnitude  of  the  emitted 
power.  The  antenna  pattern  depends  on  the  orientations, 
positions  in  space,  amplitudes  and  phasing  of  the  emitters. 
(9:109) 
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The  radiation  pattern  also  depends  on  the  distance  from 
the  antenna  that  the  field  Is  measured.  If  the  pattern  Is 
measured  sufficiently  far  from  the  antenna  for  there  to  be 
no  change  in  pattern  with  distance,  the  pattern  is  the 
far-field  pattern.  Measurements  at  lesser  distances  yield 
near-field  patterns,  which  are  a  function  of  both  angle  and 
distance.  The  pattern  may  be  expressed  in  terms  of  the 
electric  field  intensity  (field  pattern)  or  in  terms  of  the 
radiation  intensity  (power  patterns)  (6:604). 

The  pattern  usually  has  many  lobes.  The  main  lobe  is 
the  lobe  containing  the  direction  of  maximum  radiation.  Any 
lobe  other  than  the  main  lobe  is  called  a  minor  or  side 
lobe.  Typically  the  side  lobes  are  alternately  positive  and 
negative  valued.  A  pattern  in  its  most  general  form  may  be 
complex-valued.  In  that  case,  the  magnitude  of  the  electric 
field  pattern  |E(£)|  or  the  power  pattern  P(©)  is  used  to 
generate  the  antenna  pattern  (9:29). 

With  antennas  such  as  dipoles  and  horns,  the  radiation 
pattern  is  usually  expressed  in  terms  of  the  power  density 
(power  patterns).  Analyzing  array  antennas,  however, 
involves  the  addition  of  the  field  contributions  from  the 
entire  array.  The  antenna  pattern  is  then  based  on  the 
field  pattern  because  the  magnitudes  and  phases  of  the 
elements  must  be  considered  (6:627). 
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For  most  phased-array  radars  the  antenna  pattern  is 
developed  by  making  certain  geometric  constraints  such  as  a 
planar  array  with  evenly  spaced  emitters.  This  approach 
was  not  used  in  this  case  because  the  sparse  3-dimensional 
geometry  of  the  array  did  not  lead  to  such  simplifications. 
However,  the  analysis  for  the  sparse  array  was  accomplished 
in  a  similar  manner  as  with  an  n  element  array.  The  method 
is  developed  with  the  following  examples. 

The  emitters  were  assumed  to  be  isotropic  point  source 
radiators  with  the  same  amplitude  as  the  reference  source. 
In  the  first  example,  two  identical  point  source  radiators 
were  spaced  a  distance  d  apart.  In  Figure  2  the  field  was 
determined  at  a  point  p,  far  from  the  emitters  (5:398). 


Figure  2.  Two-Element  Array  of  Nondirectlonal  Emitters 


Because  the  point  p  was  very  far  from  the  emitters,  the 

distance  r^  was  approximated  by 

r  *  r  +  dcos©  ( 28 ) 

i  z 

The  sum  of  the  electric  fields  at  p  was 

E  =  E  +  E  eJv  (29) 

T  1  2 

where 

V»  =  ftdcosO  +  a 


a  =  progressive  phase  shift 
dcos©  =  path  difference 

The  progressive  phase  shift  is  used  to  electronically 
steer  the  antenna  main  beam.  For  linear  arrays,  adding  a 
constant  phase  to  each  emitter  will  point  the  beam  Into  the 
new  direction  of  interest.  For  the  two  identical  emitters 
the  broadside  team  pattern  is  obtained  when  a  =  0 . 
Therefore, 


E  =  E  1  (  1  +  ejw)  1 

(30) 

T  1  1  V  '  1 

*  E  ■/  (l  +  cosy)2  +  sln2^ 

(31) 

"  2E1cosCvy/2) 

(32) 

The  field  pattern  only  depends  on  the  distance  between 
the  two  emitters,  d,  and  the  direction  9,  When  d  =  X/2, 
then  v'  =  nc os9  and  the  field  pattern  is  shown  in  Figure  3. 

The  next  example  used  three  emitters  in  a  plane  but  not 
linearly  spaced.  Figure  4  shows  the  three  emitters  spaced 
distances  d^  and  ds  apart  (5:400). 
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The  electric  field  was  again  measured  at  point  p,  far 
from  the  array  such  that  the  distances  r^  and  r#  were 
approximated  by 

r  =  r  -  d  cos©  r  -  r  -  d  cos© 

2122  BISS 

When  emitter  1  was  the  reference  emitter,  the  total 
electric  field  was  the  sum: 

iV,  Jv' 

E«E+Ee  +Ee  (33) 

t  l  2  a 

where 

V  -  fid  cos  &  +  a 

\  l  l  l 

y/  ~  the  phase  difference 

t 


a  =  the  progressive  phase  shift 
If  all  the  emitters  are  identical  the  total  electric  field 
is  written  as : 

iv',  iv 

Et  ■  Et  (1  +  e  ♦  e  )  (34) 

The  antenna  pattern  would  be  a  graph  of  the  magnitude  of  the 
electric  field 

iv, 

Et  =  EJ  ( 1  +  e  +  e  )  I  (35) 

Because  the  positions  of  the  emitters  were  not  linear 
or  symmetrical,  the  phase  differences  could  not  be 
simplified  into  a  single  function  as  was  the  case  for  the 
two  point  emitters.  To  determine  the  magnitude  of  Et,  the 
exponent  terms  were  recast  into  real  and  Imaginary  parts. 

The  magnitude  is  given  by  the  square  root  of  the  sum  of  the 
squares  of  the  re.Jl  and  imaginary  parts. 
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*  J 

Et  =  Ef[(  1  +  cosy^  +  cosyj2  +  (  sin^  +  siny^)2]  '*  (36) 

This  method  made  use  of  the  principle  of 
superposition,  and  can  be  generalized  into  three  dimensions 
with  any  number  of  emitters ;  i  . e . ,  the  total  field  was 
the  sum  of  the  fields  due  to  each  component.  From  each 
emitter  the  phase  difference  depended  on  two  parts:  1)  the 
path  difference,  and  2)  the  phase  shift.  For  isotropic 
point  source  emitters,  the  amplitude  could  vary  so  that  in 
general,  the  total  electric  field  was  the  cum  over  all 
emitters  as  given  by  equation  (37)  (8:99-106). 

ET  =  I  EV  =  I  \  exp[  +  JvJ  (37) 

l  V 

where 

A  *=  amplitude  of  each  emitter 

=  path  difference  of  each  emitter  to  the  point  p 
*  phase  shift  added  to  point  the  beam 

The  path  difference  was  the  difference  between  the 
position  vectors  to  the  emitters  and  the  pointing  vector 
from  the  refrence  emitter,  1  .  This  difference  was 

p 

determined  by  the  projection  of  the  position  vector  to  each 
emitter  onto  the  pointing  vector,  (  r^-1^). 

The  maximum  electric  field  in  any  direction  occurs  when 
the  phase  difference  is  equal  to  zero.  From  equation  (37), 
this  occurs  when : 

/7A  «■  -yr  (38) 
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To  point  a  beam  maximum  in  a  given  direction  lp,  the  phase 
shift  is  calculated  from: 

<*(  <39> 

a 

The  direction  the  beam  was  pointed  is  defined  as  1  and 
rearranging  equation  (39)  gives: 

V  -  r  »1  1  (40) 

To  evaluate  the  electric  field  anywhere  else  with  the 
main  lobe  still  pointed  along  1  ,  the  pointing  vector,  1  , 

o  p 

only  needs  to  be  moved.  Equation  (37)  can  be  rewritten  as: 

Et  =  I  \  *XP[  ^Pfc#1p  '  VU  3  H1) 

X 

Hall-Power  Beamwldth 

A  parameter  to  describe  how  the  power  in  a  direction  is 
emitted  is  the  half-power  beamwldth  (HPBW).  It  is  the 
angular  separation  of  the  points  in  the  main  beam  where  the 
power  pattern  equals  one-half  of  the  maximum  power.  Because 
the  power  is  related  to  the  square  of  the  electric  field, 
these  points  correspond  to  the  value  1 / ' for  the  field 
pattern  |E(©)|.  The  smaller  the  HPBW  angle,  the  more 
concentrated  the  radar  energy  along  the  boresight  direction. 
(9:30) 

To  calculate  the  half-power  point  of  the  beam,  the 
magnitude  of  the  field  in  a  direction  1  was  compared  to  the 
maximum  field  strength  which  was  along  the  boresight 
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direction/  1^.  The  half-power  point  in  one  direction  from 
the  boresight  direction  was  found  when  the  ratio  of  these 
magnitudes  was  given  by: 


I  Mi,)  l 

l  Et(U  I 


0.70710678 


(42) 


The  half-power  beamwidth  angle  (HPBW)  was  determined  by 
first  assuming  the  beam  was  symmetr ica 1  about  the  boresight. 
The  pointing  vector  was  incrementally  moved  through  a  plane 
from  the  boresight.  The  electric  field  was  evaluated  at 
each  step  until  the  half-power  point  (equation  47)  was 
found.  The  HPBW  was  calculated  as  double  the  size  of  the 
angle  to  the  half-power  point. 


-Ji 


.  * 


« 
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IV.  Arrays  And.  Program 


To  analyze  the  effect  of  the  relative  motion  of  the 
satellites  on  the  beam  quality  a  computer  program  was 
developed.  The  program  used  the  locations  and 
characteristics  of  the  emitters,  determined  the  phase 
differences  and  calculated  the  magnitude  of  the  electric 
field  in  a  given  direction.  Then  by  incremental  angular 
steps  the  field  in  other  directions  was  evaluated  and  used 
to  determine  the  half-power  beamwidth.  The  HPBW  was  then 
calculated  throughout  the  hemisphere  below  the  array  each 
time  the  array  was  moved  t  seconds.  Figure  5  shows  the  flow 
chart  and  the  program  is  listed  in  the  Appendix. 

incut 

The  first  input  parameter  about  the  arrays  which  varied 
was  the  number  of  emitters  in  each  array.  The  program  was 
developed  to  put  up  to  240  satellites  in  an  array.  Arrays 
with  from  5  to  200  satellites  were  modeled  to  determine  the 
effect  of  the  density  of  emitters.  The  altitude  was  another 
variable  to  be  given  in  the  input.  Because  many  different 
operating  altitudes  had  been  possible,  the  altitude  range 
considered  was  from  500  km  to  geosynchronous  orbit. 

Other  input  parameters  were  related  to  the  emitters. 
They  were  ail  assumed  to  be  identical  point  source 
radiators.  When  calculating  the  electric  field,  since  all 
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Figure  5.  Main  Program  Flow  Chart 
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the  emitters  were  the  same,  the  amplitudes  would  factcr  out. 
For  n  emitters,  the  total  electric  field  would  be: 

Et  =  n\  £  e*p(  3^)  (43) 

i 

where 

yi  =  phase  difference  for  each  emitter 
ft  =  2 n/K 

Because  the  HPBW  is  determined  from  the  ratio  of  the 
electric  field  in  the  desired  direction  versus  the  field 
along  the  boresight  of  the  beam,  the  amplitudes  of  the 
emitters  will  cancel.  Therefore,  as  long  as  the  emitters 
are  identical,  the  amplitudes  can  be  arbitrary  and  for  the 
program  they  were  assumed  to  be  equal  to  one. 

The  wavelength  of  the  radiated  signal  is  another  input 
parameter.  The  magnitude  of  the  electric  field  is  directly 
proportional  to  the  dimensions  of  the  array  and  inversely 
proportional  to  the  wavelength  (  ri/\) .  The  units  of  r^  and 
X  therefore,  must  be  same.  In  the  program, 

M  =  3.906012  km*/s*  (44) 

For  the  mean  motion,  (n  =  */ ^ /  ■) ,  to  have  the  correct  units 
of  seconds,  the  radial  distance  r  and  the  displacements  <5r, 

O 

r  60,  and  6z  were  asaumed  to  have  the  units  of  kilometers, 
o 

To  keep  the  units  correct  and  to  normalize  the  dimensions, 
the  wavelength  was  assumed  to  be  one  kilometer.  The  actual 
dimensions  of  the  array  were  not  of  interest,  only  the 
relative  beamwidth  changes.  The  wavelength  was  therefore  a 
scale  factor  for  the  relative  positions. 
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Initial  Positions 


The  next  part  of  the  program  involved  the  input  for  the 
positions  of  the  array.  The  most  common  phased  array  radars 
are  planar  arrays  made  up  of  hundreds  to  thousands  of 
elements.  The  arrays  are  usually  fixed  to  a  building  or 
platforms  and  radiate  for  specific  areas  of  coverage. 
Locating  the  array  In  space  provides  the  opportunity  to  put 
the  array  into  a  three-dimensional  configuration.  The 
different  types  of  array  configurations  considered  for  this 
project  included  two  and  three-dimensional  arrays. 

For  the  rectangular  planar  arrays,  the  x  and  y 
coordinates  of  the  emitters  were  the  input  data.  A  scale 
factor  to  increase  or  shrink  the  overall  size  of  the  array 
and  a  tilt  angle  to  rotate  the  array  about  the  e&  axis  were 
also  given  as  input. 

The  next  type  of  planar  arrays  were  called  disks. 

These  configurations  involved  placing  elliptical  or 

spherical  rings  of  varying  numbers  of  emitters  around  the 

center.  The  shape  of  the  ring  could  be  modified  by  changing 

the  eccentricity  of  the  ring.  Cone  arrays  were  made  by 

placing  emitters  at  different  radial  altitudes. 

The  final  shapes  were  spheres  of  randomly  spaced 

emitters.  The  input  was  the  lengths  of  the  sphere  in  each 

direction.  A  random  number  generating  routine  calculated 

the  e  ,  eot  and  e  displacements.  Example  configurations 
r  v  z 

are  shown  in  Figures  6-9. 
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Figure  6.  Rectangular  Planar  Array 
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Figure  7.  19  Element,  Evenly  Spaced  Disk  Array 


26 


•I  •  •* 


Figure  8.  19  Element  Cone  Array 


Figure  9.  Randomly  Spaced  Sphere 


27 


t 


Initial  Velocities 

Once  the  initial  positions  were  given  for  any  particular 
array,  the  only  remaining  information  to  obtain  was  the 
initial  velocity  conditions.  Although  the  satellites  will 
not  remain  in  a  fixed  orientation  with  respect  to  the 
reference  satellite,  under  certain  initial  conditions  the 
motion  could  be  confined  to  a  specified  region  of  space. 

In  equations  (16)  and  (21),  the  radial  and  vertical 
displacements  have  constant  and  time  periodic  terms.  The 
tangential  displacement  equation  (24)  also  has  constant  and 
periodic  terms,  but  there  is  also  a  secular  component.  To 
keep  the  array  from  forever  spreading  apart  in  the  e^ 
direction,  the  secular  component  must  bo  made  to  equal  zero. 
From  equation  (24), 

(3r  69  +  6n<5r  1  =  0  (45) 

V  o  o  oJ 

Solving  for  the  initial  tangential  velocity: 

(  46  ) 

The  tangential  Initial  velocity  therefore  depends  on 
the  initial  radial  displacement.  By  assuming  the  only 
secular  velocity  term  equals  zero,  the  motion  of  the  array 
is  restricted  to  periodic  changes  with  the  period  of  one 
orbit.  The  array  will  spread  apart  a  maximum  distance  in 
one-half  of  an  orbit  and  then  return  to  the  initial 
positions  by  the  completion  of  one  orbit.  When  the 
conditions  in  equation  (46)  are  set,  then  the  radial 
displacement  becomes  an  oscillating  function. 


r  69  «=  -  2n6r 

O  ©  O 


J 

-1 
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I 


(47) 


6r 

6r(t)  =  {"-3<5r  -  2/  (r  60  )lcos(nt)  +  — -  sin(nt) 

Equations  (16)  and  (47)  were  used  to  develop  the 
relationship  between  the  maximum  displacements  and  the 
initial  velocity  components.  The  initial  conditions  when 
t  =  0  were: 


«5r(0)  =  6r 

and 

<5r  (  0  )  =6 r 

(48a, b) 

O 

6z(0)  =  &z 

O 

and 

6z  ( 0 )  =  6z 

O 

( 49a, b) 

•  • 

The  <5r  and  6z  are  values  to  be  determined. 

O  O 

The  maximum  displacements  occur  after  half  an  orbit  which 
is  when  y/  =  rr/2.  From  equations  (16)  and  (47),  the  maximum 
displacements  were: 

6z  -  6ro  and  62  ("/J  =  6*o  (50) 

r>  r» 

The  initial  velocities  were  then  found  in  terms  of  a  given 
maximum  allowable  displacement  in  each  direction. 

•  • 

6r  =  n6r  6z  =  n6z  ( 51 ) 

o  max  o  max 

In  the  program  61  and  <5z  were  set  equal  to  one  for  the 

max  max 

bas ic  analysis . 

fieamwldth  Subroutine 

The  subroutine  EFIELD  was  used  to  evaluate  the  electric 
field  by  incremental  steps  in  the  hemisphere  below  the 
array.  The  magnitudes  of  the  electric  fields  were  used  to 
determine  the  HPBW  angle.  The  routine  used  the  locations  of 
the  emitters  to  calculate  the  phase  difference  for  each 
direction.  The  coordinate  system  used  to  determine  the  look 
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directions  and  phase  differences  is  shown  in  Figure  10.  The 
angle  4>  ranged  from  0*  to  90*  while  the  angle  ©  ranged  from 
0*  to  360*.  The  angular  steps  were  15*  in  both  directions. 
Fifteen  degrees  was  chosen  because  it  provided  enough  change 
when  calculating  the  field  without  creating  too  much  data 
for  the  plotting  program. 
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SIDE  TOP 


Figure  10.  Pointing  Vector  Coordinate  System 
To  determine  the  half-power  point,  the  point  direction, 
1  was  changed  in  the  <p  direction.  The  field  in  that 

p 

direction  was  compared  to  the  beam  maximum  field  value  until 
the  half-power  point  was  found.  The  algorithm  used  a 
step-down  process  to  locate  the  bounds  for  a  bisection 
search  which  determined  the  half-power  point.  The  beam  was 
assumed  to  be  symmetrical  about  the  boresight  so  that  the 
half-power  beamwldth  was  twice  the  angle  to  the  half-power 
point.  The  flow  chart  for  the  subroutine  EFIELD  is  shown 
in  Figure  11. 
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QxhlL  SubroviUng 

Once  the  half-power  beamwidths  in  all  directions  were 
calculated  and  recorded,  the  next  step  was  to  move  the  array 
forward  in  time  with  the  subroutine  MOVEIT.  The  amount  of 
minutes  to  move  was  determined  in  the  main  program  and  then 
converted  to  seconds  in  the  subroutine.  Time  was  made 
<J  linens  ionless  by  y  =  n  x  t.  The  matrices,  equations 
(27a-d),  were  evaluated  and  used  to  calculate  6IM  t ) 
and  6V  (t).  The  new  positions  <5R  ( t )  were  then  returned  to 

i  i 

the  main  program  for  the  half-power  beamwidths  to  be 
recalculated . 
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y_i.  Result  a  And.  Analysis 


Program  CIumlL 

To  test  the  validity  of  the  program,  a  linear  array  of 
20  emitters  was  used  with  a  spacing  of  X./2.  The  values  for 
the  HPBW  were  computed  and  compared  to  the  known  solutions 

A 

In  reference  6.  The  HPBWs  In  the  plane  are  shown  In 
Figure  12.  The  pattern  Is  shown  only  In  the  plane  of  the 
array  because  the  field  of  the  linear  array  Is  symmetric 
about  the  axis  of  the  array.  The  actual  three  dimensional 
antenna  pattern  Is  a  disk  shaped  figure  obtained  by  rotating 

a 

the  pattern  about  the  e0  axis  (6:63b). 

The  calculated  broadside  and  endflre  HPBWs  were  5.08* 
and  34.25*.  The  text  values  were  5.1*  and  34*.  As  Figure 
12  shows,  the  HPBW  Increased  as  the  beam  was  pointed  towards 
the  endflre  direction  (6:635). 

To  check  the  motion  subroutine,  a  random  sphere  with 
three  elements  was  used.  The  initial  positions  were 
generated  by  the  program.  With  these  positions  and  6x 

max 

and  6z  equal  to  one,  the  initial  velocities,  n,  and 

man 

the  period  were  calculated.  The  positions  at  different 
points  In  the  orbit  were  calculated  with  a  hand  calculator. 
These  values  and  the  program  results  for  one  point  are 
listed  in  Table  I.  The  program  gave  the  same  results. 
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Method 

6x 

r  66 

o 

6z 

Program 

0 

11  ^  3 

2  n 

-91.7001 

46.7161 

-91.7001 

37.1940 

193.0233 

37  .  19  41 

31.2645 

-14.7662 

31.2645 

Calculat . 

0 

-91.7001 

37 . 1940 

31.2645 

i^ErllrlLir^K. 

192 . 1637 

-14  .7662 

3 

2n 

-91 . 7001 

37.1933 

31.2645 
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Initial  Results 


Overview  Each  array  was  evaluated  at  Its  initial 
positions  and  then  after  each  time  it  was  moved.  Before  and 
after  configuration  and  HPBW  plots  were  made  for  each  array. 
The  dimensions  on  both  axes  of  the  configuration  plots  are 
In  terms  of  wavelengths.  The  figures  show  the  relative 
displacements  between  the  emitters.  All  displacements  were 
normalized  by  setting  the  wavelengths  equal  to  one.  The 
configuration  figures  show  top  views  of  the  arrays. 

After  the  half-power  beamwidth  was  calculated,  the 
angles  <f>  and  0  were  converted  to  rectangular  x  and  y 

A  A 

coordinates  in  the  e_  -  e  plane.  The  positive  z  axis  value 

z 

was  the  HPBW  value  from  the  4>,  0  direction  below  the  array. 
To  graph  the  data,  the  file  of  x,  y,  z  data  was  then  used 
with  the  SURFER  graphics  package  to  generate  contour  plots. 
When  making  the  plots,  most  default  settings  In  the  SURFER 
programs  were  used. 

In  the  initial  analysis,  the  arrays  were  only  moved  1/3 
of  an  orbit.  The  one  third  orbit  point  was  chosen  to  avoid 
the  possibility  of  masking  changes  in  the  HPBW  caused  by 
symmetric  motion  that  might  have  occured  at  the  1/4  or  1/2 
orbit  points.  The  1/3  orbit  point  was  also  chosen  to  give 
enough  time  for  a  significant  change  to  occur. 

Four  general  types  of  arrays  were  developed  for  the 
initial  analysis.  The  arrays  were  1)  a  rectangular  plane, 

2)  a  disk,  3)  a  cone  and  4)  a  sphere.  These  arrays  were 
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evaluated  to  determine  the  most  feasible  configuration  for  a 
distributed  array.  One  type  of  array  was  chosen  and  then 
modified  to  determine  the  effects  of  changing  the  initial 
conditions . 

Planar  Array  A  twenty-five  element,  evenly  spaced, 
planar  array  was  the  first  configuration.  The  top  view  of 
the  relative  positions  before  and  after  moving  1/3  an  orbit 
are  shown  in  Figure  13.  Figure  14  shows  the  HPBV  contours 
before  and  after  1/3  an  orbit.  As  with  the  linear  array, 
the  HPBW  is  good  near  broadside  and  degrades  towards  the 
sides  of  the  array.  Planar  arrays  have  limited  usable 
directions  because  the  beam  width  increases  when  pointing 
the  beam  in  the  endfire  direction. 

Disk  Array  Another  form  of  a  planar  array  was  a  19 

element,  slightly  elliptical,  disk  array.  The  emitters 

were  evenly  spaced  in  rings  about  the  center.  Figure  15 

shows  the  relative  positions  of  the  configuration.  The 

motion  of  this  and  the  rectangular  plane  are  similar  because 

they  were  both  initially  in  a  flat  plane  with  no  er 

displacements.  As  these  arrays  orbit,  the  dominant 

oscillating  motion  was  in  the  e  direction.  The  motion  in 

z 

the  e&  direction  was  considerably  smaller  than  the  initial 
displacements  and  cannot  be  seen  on  these  figures  because  of 
the  scale.  Figure  16  shows  the  HPBW  contours.  The  results 
are  similar  to  the  rectangular  planar  array  In  that  the  beam 
width  Increases  towards  the  endfire  conditions. 
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16.  Disk  Array  HPBW  Contours 


Cone  Array  The  cone  array  was  made  by  putting  the 
rings  o£  the  19  element  disk  array  at  increasing  radial 

levels.  Figure  17  shows  the  relative  positions  o£  the  cone 
array.  Because  of  the  displacements,  the  elements  moved 
more  in  the  direction  than  with  the  disk  array. 

Figure  18  shows  the  HPBW  contours.  As  with  the  planar 
arrays,  the  KPBW  increased  when  the  beam  was  pointed  away 
from  the  downward  radial  direction. 

Sphere  Instead  of  a  specific  configuration,  a  purely 
random  displacement  was  tried.  A  sphere  of  twenty  emitters 
was  used.  The  random  sphere  was  chosen  because  it  would 
reflect  the  most  general  orbit  placement.  The  satellites 
could  be  released  from  a  booster  in  the  general  area  and 
given  enough  av  to  ensure  adequate  spacing.  The  positions 
of  the  emitters  in  the  sphere  are  shown  in  Figure  19.  The 
HPBW  contour  plots  are  shown  in  Figure  20. 

Figure  20  shows  the  sphere  had  very  similar  HPBWs  in 
all  directions  below  the  array.  Unlike  the  planar  arrays, 
the  sphere  would  not  be  limited  to  certain  angles  of 
operation.  For  a  three-dimensional  array  in  space,  the 
capability  to  have  the  same  beamwidth  in  all  directions 
would  be  ideal.  The  array  could  be  used  to  track  targets 
almost  from  horizon  to  horizon. 
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Figure  17.  Cone  Array  Configurations 


Figure  18.  Cone  Array  HPBW  Contours 
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Sphere  Analysis 


Because  the  sphere  had  good  beamvidths  in  almost  all 
directions,  it  was  chosen  to  use  for  further  analysis.  The 
next  step  was  to  modify  different  spheres  to  determine  the 
effects  of  changing  the  various  input  parameters.  Although 
there  were  many  input  parameters,  only  a  few  actually  caused 
the  HPBW  to  change.  The  input  parameters  that  were  varied 
were  the  altitude,  wavelengths,  number  of  emitters  and 
relative  positions  of  the  emitters. 

The  altitude  was  found  to  have  no  effect  on  the  HPBW. 
One  array  was  run  with  varying  altitudes  and  moved  1/3  ar 
orbit  each  time.  The  relative  positions  and  HPBW  contours 
were  the  same  in  each  case.  This  result  was  correct 
because  in  equation  (41),  r  was  not  a  factor  in 

O 

determining  the  electric  field.  Also,  equations  (27)  show 
that  the  altitude  was  not  a  variable  and,  therefore,  did  not 
effect  the  relative  positions.  The  altitude  was  therefore, 
eliminated  as  a  limiting  parameter. 

The  wavelength  had  a  direct  effect  on  the  HPBW.  The 
half-power  beamwidth  could  be  approximated  by  equation  (52) 
because  the  radar  analysis  was  in  the  far-field  and  the 
configurations  were  sparsely  spaced  arrays  with  tne  aperture 
dimension  D,  much  greater  than  one  wavelength.  Equation 
(52)  yields  the  HPBW  in  radians  (7:362-363). 

HPBW  %  X/D  (52) 
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Equation  (52)  s*~ovs  the  HPBW  was  inversely  proportional 
to  the  size  of  the  array.  For  a  given  wavelength,  the 
larger  the  array,  the  set  ' .  r  the  HPBW.  The  relative 
dimensions  were  normalized  to  wavelengths  by  setting  the 
wavelength  equal  to  one. 

Two  spheres  with  different  radial  lengths  were  run  to 
show  the  effect  due  to  the  size  of  the  array.  The  initial 
radius  of  sphere  1  was  75X  and  the  radius  of  sphere  2  was 
200\.  The  initial  and  1/3  orbit  cases  were  run  and  Table  II 
lists  the  minimum  and  maximum  HPBW  values.  Figures  21  and 
22  show  the  HPBW  contours.  As  expected,  the  size  of  the 
HPBW  decreased  with  the  larger  array.  To  show  the 
similarity,  the  contour  plots  for  sphere  2  were  scaled  to 
twice  the  values  of  sphere  1.  Although  the  values  of  the 
HPBWs  changed,  the  shapes  of  the  contours  did  not. 


Table  II.  Minimum  and  Maximum  HPBW  Values 


Sphere 

Half-power  Beamwidth  (deg) 

Minimum 

Maximum 

1  (0) 

0 . 2724 

0.4570 

1  (1/3) 

0.1639 

1.6141 

2  (0) 

0.1022 

0.1714 

0.0615 

0.6005 
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Figure  21.  Sphere  1  HPBW  Contours 
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MOVED  1/3  ORBIT 


Figure  22.  Sphere  2  HPBW  Contours 


49 


For  all  the  spheres  there  were  spikes  in  the  HPBW 
contour  plots  of  the  1/3  orbit  configurations.  A  possible 
cause  for  these  spikes  was  the  specific  placement  of  the 
emitters.  Another  cause  might  have  been  the  number  of 
emitters  in  the  sphere.  A  third  possibility  was  the  spikes 
were  just  the  result  of  the  orbit  motion. 

For  the  first  case,  three  other  spheres  with  the  same 
number  of  emitters  were  run  but  with  different  random  number 
seeds.  The  seed  number  caused  the  random  generating  routine 
to  make  different  initial  positions. 

The  HPBW  contour  plots  for  spheres  A,  B,  and  C  are  in 
Figures  23,  24,  and  25.  The  scaling  on  the  contour  maps  was 
not  all  the  same.  Some  plots  were  scaled  as  much  as  20 
times  to  enhance  the  small  changes  in  the  beamwidth.  The 
actual  minimum  and  maximum  HPBW  are  given  in  Table  III. 


Table  III.  Three  Random  Sphere  Beamwidths 


Sphere 

Half-power  Beamwidth  (deg) 

Before 

Minimum 

Maximum 

A 

0.2863 

0.3353 

B 

0 . 2092 

0.4841 

C 

0 .2337 

0.4869 

After 

A 

0.1458 

1 . 3013 

B 

0.1112 

1.5616 

C 

0 . 2079 

1.0755 
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Figure  23.  Random  Sphere  A  HPBW  Contours 
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Figure  24.  Random  Sphere  B  HPBW  Contours 
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Figure  25.  Random  Sphere  C  HPBW  Contours 
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The  HPBW  figures  show  each  array  was  initially  unique. 
Spheres  A  and  B  had  larger  HPBW  values  in  the  forward  e0 
direction.  After  moving  1/3  an  orbit,  the  beamwldth 
improved  in  the  forward  direction  and  was  worse  towards  the 
rear  edges.  All  the  spikes  appeared  in  the  general  area  to 
the  left  rear  of  the  center. 

The  final  modification  was  to  determine  the  effect  of 
the  density  of  the  sphere  on  the  HPBW.  Three  other  spheres 
with  different  numbers  of  emitters  were  run.  The  spheres 
had  5,  25,  and  200  emitters.  Figures  26,  27,  and  28  show 
the  HPBW  contours.  The  figures  show  the  beamwidth  depended 
strongly  on  the  number  of  emitters  in  the  sphere.  In 
certain  directions,  the  HPBW  increased  because  of  the 
positions  of  the  emitters.  In  those  directions,  the 
projected  aperture  dlmenslcn  D  decreased. 

As  Figure  26  shows,  the  5  element  sphere  was  poor  in 
many  directions.  When  the  number  of  emitters  was  increased, 
the  configuration  more  resembled  a  solid  sphere  with  the 
projections  in  all  directions  being  nearly  Identical.  The 
HPBWs ,  therefore,  were  small  and  similar  in  all  directions. 
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Figure  27.  Twenty-five  Element  Sphere  HPBW  Contours 


gure  28.  200  Element  Sphere  HPBW  Contours 
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For  each  input  modification  analyzed,  the  spikes  were 
present.  The  spikes,  therefore,  appeared  to  have  been 
related  to  the  motion  of  the  array.  Another  sphere,  was 
used  to  show  the  change  of  the  HPBW  at  several  points  in  the 
orbit.  The  sphere.  Sphere  X,  nad  20  elements.  The  HPBW 
contour  plots  were  generated  for  each  l/6in  of  an  orbit 
and  are  shown  in  Figures  29  -  32. 

The  figures  show  the  degraded  HPBW  began  in  the  forward 

A  A 

eQ  direction  and  progressed  to  the  rear  eQ  direction  as  the 
array  moved  through  one  half  of  an  orbit.  After  half  of  an 
orbit  the  HPBW  values  returned  to  the  original  0.24  -  0.27 
degree  range.  However,  the  humps  moved  indicating  the 
directions  of  poor  beamwidths  made  a  180  degree  shift  from 
the  front  to  the  rear.  The  reason  tor  tnis  shift  was  that 
after  half  an  orbit,  the  entire  array  had  reversed  itself. 

A?  the  array  returned  to  its  original  positions,  the  HPBW 
values  repeated  the  patterns  in  moving  from  the  front  to  the 
rear  . 
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Figure  29.  Early  Orbit  HPBW  Contours  for  Sphere  X 
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MOVED  2/3  ORBIT 
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Figure  31.  Late  Orbit  HPBW  Contours  for  Sphere  X 
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Figure  3 
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Final  HPBW  Contour  for  Sphere  X 
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Summar.y 

The  most  useful  conf igura t i on  was  found  to  be  a  sphere 
of  randomly  spaced  emitters.  The  sphere  provided  useful 
half-power  beamwidths  in  all  directions  below  the  array. 

The  planar  arrays'  regions  of  useful  beamwidths  were  limited 
to  just  below  the  configurations.  Because  the  random  sphere 
could  essentially  look  from  horizon  to  horizon/  it  would 
require  fewer  emitters  to  cover  a  large  area  than  a  planar 
array.  The  random  sphere  of  emitters  would  be  easier  to 
construct  than  the  other  arrays  because  the  initial 
positions  would  be  less  stringent. 

The  value  of  the  HPBW  in  a  given  direction  depended  on 
th  lumber  and  relative  placement  of  the  emitters.  The 
n^i.t  r  of  necessary  emitters  was  not  large.  The  expected 
value  of  beamwidth,  A/D,  was  achieved  in  all  directions  with 
about  twenty  emitters.  Five  emitters  was  too  few  and  200 
was  overkill. 

An  array  was  assumed  to  be  in  two-body  motion.  To  keep 
the  emittors  from  forever  drifting  appart,  there  was  one 
limiting  assumption.  The  initial  velocity  in  the  eQ 
direction  had  to  be  proportional  to  the  initial  radial 
displacement  of  the  emitter,  (equation  46).  As  an  array 
progressed  in  its  orbit,  the  oscillatory  motion  caused  the 
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directions  of  poor  beamwidths  to  predictably  shift.  The 
half-power  beamwidths  actually  improved  in  certain 
directions  and  only  slightly  degraded  in  others. 

Recomme  nd  a  1 1  ons 

Although  only  spheres  were  analyzed  for  this  thesis, 
other  three-dimensional  shapes  could  be  evaluated.  The 
regions  of  antenna  coverage  could  be  varied  by  using 
ellipsoids  with  different  axis  lengths.  Certain 
orientations  of  the  ellipsoid  could  improve  the  half-power 
beamwidths  but  the  effect  of  the  motion  must  be  considered. 

For  this  thesis,  the  emitters  were  assumed  to  be  point 
source  radiators.  Follow  on  work  could  use  other  types  of 
space  qualified  radar  emitters  to  predict  the  effects  of 
orbital  motion  on  the  capabilities  of  the  arrays. 
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Computer  Pr  oqram 
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c 

c 

25 

c 


prog  rani  newpat 

integer  1 1 1 ,  w 

real  nn,mu,  min,  lam,  Eo 

character*12  Outfile,  infile 

common  ro,n, lam, nn, mu, radius, alt, Eo 

dimension  R(240,3),V(240,?),ro(240,3),vo(240,3) 

wr ite ( *, 5  ) 

f ormat  ( lx,  ’ ENTER  THE  NAME  OF  THE  LOCATION  OUTPUT  FILE') 
read  ( *, 10  )  Outf  i  le 
format  (al2) 

open  (unit-2, file=Outfile, status-' new'  ) 
wr  i  te  (  * ,  1 2  ) 

format (lx, 'ENTER  THE  NAME  OF  THE  FILE  FOR  INPUT’) 
read  (*,10)  infile 

open  (unlt-5,file-infilfe,status= 'old' ) 

THE  NUMEER  OF  ARRAY  ELEMENTS  IS  N 
read (5,15)  n 
wr  i  te  (  2, 1  5  )  ri 
format  ( i  4  ) 

INPUT  THE  WAVELENGTH  (lam)  IN  KILOMETERS 

read  (5,20)  lam 
format(fl5.7) 

INPUT  THE  ALTITUDE  OF  THE  CENTER  IN  KILOMETERS 
read (5,20)  alt 

INPUT  MAXIMUM  DISPLACEMENTS  IN  Er  AND  Ez  DIRECTIONS 
THESE  MUST  BE  IN  KILOMETERS 
read (5,20  Jdrmax 
read  (  5, 2C  Idzmax 

radius  -  6378.145  +  alt 
mu  =  3. 98601 2e 05 
vcir  =  sqr t ( mu/rad i  us  ) 
nn  =  sqr t ( mu/ ( radi us**3 ) ) 
period  =  2  .  *3 . 14 1 59  26 5 3 5/nn 

SPECIFY  THE  TYPE  OF  ARRAY  TO  USE' 

PICK  ONE  OF  THE  FOLLOWING  CASES' 

1)  RECTANGULAR  PLANE' 

2)  HEXOGONAL  PLANE’ 

3)  CIRCULAR  PLANE' 

4)  ELLIPTICAL  PLANE' 

5)  RANDOMLY  SPACED  ELLIPSOID' 
read  (5,25)  i ck 

f  or  ma  t  ( i  4  ) 
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DETERMINE  THE  SPACING  OF  THE  ELEMENTS  WR2  THE  CENTER 


if  (ick. 

eq . 1 )  then 

ca  1 1 

rectangle 

else  if 

(  ick . eq . 2 ) 

then 

call 

hexogon 

else  if 

(  ick . eq . 3 ) 

then 

call 

c  i  r  c 

else  if 

(  ick . eq . 4 ) 

then 

call 

platter 

else  if 

(  ick . eq . 5 ) 

then 

cal  1 

blob 

end  if 
c 

min  =  0.0 
wi i te ( 2 , 3 5  )  min 
35  £ormat(el5.7) 

do  50  i  =  1 ,  n 

write(2,40)i,ro(i,l),ro(i,2),ro(i,3) 

40  £ormat(lx,i4,3(2x,el5.7)) 

50  continue 

c 

wr i te ( * , 55) 

55  format ( I x, ' INITIALIZING  R()  FROM  RO()') 

do  65  i  =  1 ,  n 

do  60  j  =  l,  3 
r  (  i ,  j  )  =  r  o  (  i ,  j  ) 

60  continue 

65  continue 

c 

wr  i  te  (  *  ,  7  0  ) 

70  format ( lx, ’ INITIALIZING  VO ( ) ' ) 

do  7  5  i  =  1 ,  n 

vo(i,2)  =  -2 . 0*nn*r o ( i , 1  ) 

7  5  cont inue 

do  80  i  =  l,n 

vo  (  i  ,  3  )  =  nn*dzmax 

vo(i,l)  =  nn*(drmax  -  4  .  0*r o ( i , 1 )  ) - ( 2 . 0*vo ( i , 2 ) ) 
80  continue 

c  write(2,85)  alt 

85  format ( lx, ' ALTITUDE  =  ',el5.7) 

c  wr ite ( 2, 90 )  drmax 

9  0  for  mat ( lx, ' DRMAX  =  '  ,e!5.7) 

c  write(2,95)  dzmax 

95  format ( lx, ' DZMAX  =  ’,el5.7) 

c  do  97  i-l,n 

c  write(2,96)  i , vo ( i , 1 ) , vo ( i , 2 ) , vo ( l , 3 ) 

c9  6  for  ma  t  (  lx,  i  4  ,  3  (  2x,  e  1  5 . 7  )  ) 

c97  continue 

delmin  =  period/3.0/60.0 
do  125  boo  -1,2 
call  efield(r,min) 
min  =  min  +  delmin 
call  move  it ( r , v, min, vo ) 
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100 

c 

c  10  5 


107 
110 
12  3 


1  4  A 
145 


c 


200 

22C 

c 

c 

c 

c 

2  30 


c 

c 

c 

c 

c 


240 


c 


wr  ite ( 2, 100 )  min 
format(el5.7) 
write(2,105) 

format ( lx, 'NEW  POSITIONS'/) 
do  110  i  =  l,n 

write(2<107)i,r(i,l),r(i,2),r(i,3) 
format(lx, i4,3(2x,el5.7)  ) 
cont i nue 
continue 
post  =  9999.0 
wr ite ( 2, 144  )  post 
wr i te ( * , 1 4  5 } 
format(el5.7) 
i ormat ( lx,  'FINISHED  '  ) 
close ( unit=S ) 
close (unit =2 ) 
stop 
end 

subroutine  rectangle 
common  ro,n,lam,nii,mu,radius,alt,Eo 
dimension  a(240),b(240),ro(240,3) 
real  lam 

do  220  i = 1 , n 
read (5,200)  a ( i ) , b( i ) 
format  (  2  f  1 0 . 5  ) 
continue 


INPUT  THE  TILT  ANGLE  FOR  THE  PLANE  (DEGREES)' 

TOP  GOING  AWAY  IS  POSITIVE’ 

read  (5,230)  alpha 
format(flO.S) 

alpha  =  alpha* 3 . 1  4 1  59/1 8 0 . 0 
ss  =  s i n ( a lpha  ) 
cc  -  cos(alpha) 

DISTANCES  ARE  MULTIPLES  OF  THE  SIGNAL  WAVELENGTH 


ENTER  SCALE  FACTOR  (d)  FOR  SPACING  (2*LAM,  ETC)’ 
read  (5,230)  scale 


do  240  i  =  1 ,  n 

ro(i,l)  =  b ( i ) *ss*scale*lam 
ro(i,2)  =  -a ( i ) *sca le *  lam 
ro(i,3)  =  b( i ) *cc*scalfc*lam 
cont inua 
return 
end 


subroutine  hexogon 

return 

end 

subroutine  circ 
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return 

end 

c 

subroutine  platter 

common  ro,n, lam, nn, mu, radius, alt,Eo 
dimension  10(240,3) 
integer  nrings 
psi  =  0.0 
j  =  0 

c  A  AND  DEPTH  MUST  BE  IN  KILOMETERS 

read(5,250)  a , ecc , n r  i ngs 
250  format(2(fl0.3),i5) 

do  260  i  =1, nr  ings 
read(5,252)  r i ng fac , i nr , dept h 
252  £ormat(fl0.3, i 5 , f 1 0 . 3) 

dpsi  =360.0/ inr 
ao  -  a  *  ring fac 
tio  =  au  *  sqrt(l.-ecc**2) 
do  265  d  =  1 , i nr 
j  =  j  + 1 

ppsi  =  psi  *  3.14159/180.0 
r  o  (  }  ,  2  )  =  -ao*cos ( pps i ) 
ro(j,3)  =  bo*s in ( ppsi ) 
r  o ( j , 1 )  =  depth 
psi  =  psi  +  dpsi 
265  continue 

260  continue 

return 
end 
c 

subroutine  blob 

common  r  o, n,  lam, nn , mu, rad i us , a  1 1 , Eo 
dimension  ro(240,3) 

real  legr,legthe,legr,numhl,numb2,numb3/jj,kk,jo/ko 
integer  y , w 

C  THIS  ROUTINE  GENERATES  RANDOM  LOCATIONS 

c  FOR  THE  ARRAY  ELEMENTS 

c 

c  INPUT  MAXIMUM  LENGTHS  OF  EACH  AXIS  OF  THE  ELLIPSOID 

c 

c  THE  AXIS  LENGTHS  ARE  IN  KILOMETERS" 

c 

C  Y  IS  ANY  ODD  NUMBER  BETWEEN  0  AND  67108863 

c 

read  (5,300)  y 
300  £ormat(i8) 

read  ( 5, 311 )  legr 
read  (5,311)  legthe 
read  (5, 311 )  legz 
311  format  (  f  15.  7  ) 
c  establishing  the  center 
r o  ( 1 , 1 )  =  0.0 

ro ( 1, 2 )  =  0.0 
rod, 3)  =  0.0 
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do  320  i  -  2, n 

numbl  -  rand(y) 

ro(l,l)  =  legr  *  (2*riumbl~l) 

numb 2  =  rand(y) 

ro(i,2)  =  legthe  *  (2*numb2-l) 

numb3  =  rand(y) 

ro(i,3)  =  legz  *  (2*numb3-I) 

320  continue 

return 
end 
c 

function  rand(iy) 
c 

c  pseudo  random  number  generator  on  interval  (0,1) 

c  Collected  Algorithms  of  the  CCM  #266 

c  assumes  2**31  integer  math 

c  iy  is  odd  integer  between  0  and  67106863 

c  don't  change  it  after  first  call 

c 

dimension  k(3) 
data  k/ 25, 25,5/ 
c 

do  400  i  =  1,3 
iy  =  k ( i  )  *  iy 

iy  =  iy  -  ( iy/67108864  )  *67108664 
400  continue 

rand  —  r  e  a  1 ( iy)/671088G4.0 
return 
end 
c 

subroutine  e f i e Id ( r , mi n ) 
common  r o , n , lam, nn , mu , r ad i us , a  1 1 , Eo 
real  lpr,lpz,lor,loz/lot/lpf,jj,kk,jo/ko/lam,mi:. 
real  aa,bb/zz,diffp,diffa,ea,ep/chek,mu,nn 
character*12  out  name, out  file 
dimension  doto(240),dotp(24C),direct(40,80) 
dimension  r o ( 240 ,  3 ) , r ( 240, 3 ) 
c  wr  i  te  (  *  ,  500  ) 

c600  for  mat ( lx, ' Enter  the  name  of  the  Efield  output  file') 

c  read (*,510)  out  name 

c510  £ormat(al2) 

c  open  (unit- 7,  f i le= ' f ield . out ' , status= ' new' ) 

open  (unit=8,file='plotfile',status='new' ) 
c  wr  i  te  (  7 , 51 5  )  .iin 

c  write(8,515)  min 

515  format.  (Ix,el5.7) 
c 

C  EACH  ELEMENT  IS  ASSUMED  TO  BE  AN  ISOTROPIC  POINT 

c  SOURCE  EMITTING  ITS  ENERGY  IN  PHASE 

c 

c  ALL  ELEMENTS  EMIT  THE  SAME  IN-PHASE  AMPLITUDE,  Eo 

Eo  =  1.0 
i count  =  0 
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itag  =  1 

do  600  3=0,90,15 
j  j  -  3 
jo  =  j  j 

c  this  needs  to  go  0  to  360 
do  5S0  k  =  0 ,  360,  30 
suml  =  0.0 
s  um2  =  0.0 
kk  =  k 
ko  =  k  k 

i£(k.eq.360)  then 
goto  590 
end  if 

phi  =  jj  *  3 .1415926535/1 80 . 0 
the  =  kk  *  3.1415926535/180.0 
lor  =  -sin(phi) 
lot  =  cos ( ph i ) *cos { the } 
loz  -  cos ( ph 1 ] *s  i  n ( the  ) 
lpr  =  lor 
lpt  =  lot 
lpz  =  loz 
do  520  i -1, n 

doto(i)  =  (r(i,l)*lor  +  r(i,2)*lot  + 

r(i,3)*loz)*6. 2831853/lam 

dotp(i)  -  (r(i,l)*lpr  +  r(i,2)*lpt  + 

r  < i , 3) *lpz ) *6 . 2831853/lam 

suml  =  suml  •»  cos  ( dotp  ( i  ) -doto  ( i  )  ) 
sum2  =  surn2  +  s  i  n  ( do t p  {  i  ) -dot  o  (  A  )  ) 

520  continue 

emax  =  Eo *sqr t ( s uml *s uml  +  sum2*sum2) 

S  1  -  0.0 
sum 2  =  0.0 

c  *******  STEP  DOWN  SEARCH  ******** 

delta  =  1.0 
jcount  -  0 

999  i £( jcount . ge . 360 )  then 
band  =  1.0  -  ratio 
i f ( Land . 1 t . 0 . 000000001 )  then 
hpbwl  -  360.0 
else 

hpbwl  =  5555.5555 
end  if 
j  j  =  jo 
jcount  =  0 
goto  1001 
end  if 
sum2  =  0.0 
suml  =  0.0 

phi  =  jj  *  3.1415926535/180.0 
the  =  kk  *  3.1415926535/180.0 
lpr  =  -sin (phi' 
lpt  =  cos ( ph i ) *cos ( the ) 
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1  p  z  =  cos ( ph i ) *s i n ( t he  ) 
do  530  i  =  1,11 

dotp(i)  -  (r(i,l)*lpr  ■»  r(i,2J*lpt  + 

r(i,3)*lpz)*6. 2831853/lam 

suml  -  suml  +  cos  (  d  o  tp  (  i  )  -dot. u  (  i  )  ) 
sum2  =  suni2  +  s  i  n  (  d  o  t  p  (  i  ) -dot  o  (  i  ) 

530  continue 

ENOW  =  Eo*sqr  t  (  suml *suml  +  sum2*sum2) 
suml =0.0 
s  um2  =  0.0 

ratio  =  ENOW/emax 
differ  =  ratio  -  0.7071 
i f ( d i f f er . gt . 0 . 0  )  then 
j  j  =  jj  +  delta 
jcc^nt  =  j count  +  1 
goto  999 
end  if 

c  END  STEP  DOWN  SEARCH 
c 

c*********************************** 

C  BISECTION  SEARCH  FOR  HPBW1  BETWEEN  JJ-1.0  AND  JJ 

c  DOING  THIS  FOR  JJ  ONLY  NOW 

1155  aa  =  j  j  -  delta 
bb  =  j  j 
tul  =  0.00001 
no max  =  Li 
icount  -  1 

3333  i f ( icount . It . nomax )  then 
flag  =  1 

2  2  =  a  a 
goto  222 

5555  EA  =  ENEW 
7777  p  =  aa  +  (bb  -  aa)/ 2.0 
,.z  =  p 
flag  =  2 
goto  222 
4444  EP  =  ENEW 

diffp  =  ((EP/emax)  -  0.70710676) 
chek  =  abs ( bb-aa  ) 

if(abs (diffp) . It . t  o 1 .or.chek.lt .tol )  then 
hpbwl  =  2 . 0*abs ( jo-p ) 
j  j  -  jo 
icount  =  0 
goto  1001 
e  nd  i  f 

diffa  =  ( (EA/emax)  -  0.70710678) 
if(diffa*diffp.gt.0.0)  then 
aa  -  p 
else 

bb  =  p 
end  if 

icount  =  icount  +  1 
goto  3333 
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else 

goto  8088 
end  if 

£  *********************************A*** 

222  suml  -  0.0 
sum2  -  0.0 

phi  =  zz  *  3 . 1 41 0928030/180 . 0 

the  -  kk  *  3.1415926S35/180.0 

1 pr  =  -s i n ( ph i  ) 

lpt  =  cos ( ph i ) *cos ( the  ) 

ipz  -  cos ( ph I ) * s i n ( the ) 

do  5  50  i  -  1,11 

dotp(i)  =  (r(i,l)*lpr  +  r(i,2)*lpt  + 

r  [  i  ,  3  )  *lpz )  *6  . 2831853/ lam 
suml  =  suml  +  cos ( dotp ( i ' -Joto ( i 1 ) 
suml  -  suml  +  s i n ( dotp ( i ) -doto ( i ) ) 

550  continue 

ENEW  -  Eo*sqr t ( suml**2  +  sum2**2) 

r-  ******************************* 

if (flag. eg. 1)  then 
goto  5555 
end  if 

i f ( f  lag . eg . 2 )  then 
goto  ‘3444 
end  if 

c  DIDN'T  WORK 

8886  hpbwl  -  2222.2222 

£  ****************************** 

1001  continue 

c  write(7,560)  jo,ku,hpbwl 

560  format(lx,2(2x,f6.1),2x/el5.7) 

c  sending  data  to  the  file  to  make  the  plot 
tto  =  ko*3. 1415926535/ 180.0 
ppo  =  jo*3 . 141  59  26535/180 . 0 
XXX  1  cos ( ppo ) *cos ( t to ) 

yyy  -  cos ( ppu ) i n ( t t o ) 
zzz  =  hpbwl 

write(8,57G)  xxx,yyy,zzz 
570  format(l>./3(2x/el5.7)) 

D  j  =  jo 

590  continue 

c  wr  i  te  (  *  ,  59  5  ) 

c  59  5  format ( lx, ’WORKING’  ) 

600  cont inue 

wr i te ( * , 610  ) 

610  format ( lx, 'FINISHED  DOING  THE  E-FIELD') 
c  close  ( un i t  =  7  ) 

close  ( unit  =  8  ) 
return 
end 
c 
c 

subroutine  move i t ( r , v, mi n, vo ) 
common  ro,n,  lam,nn,mu,radius,alt/Eo 
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dime  ns  ion  phirr ( 3, 3 ) ,  ph i r  v ( 3, 3 ) , ph i vr ( 3, 3 ) , phi vv ( 3, 3 

dimension  to(240/3),r(240#3)/vo(240#3)/v(240/3) 
real  lam,  rm,  mu,  mi  n 
c 

c  wr i te ( * ,  700  ) 

c  7  0  0  formal.  (  lx,  ’  INPUT  THE  TIME  IN  MINUTES  TO  MOVE’) 

c  r  ead i * , 710 )  min 

c710  formaU  f  IS.  7  ) 

t  =  6  C  .  0  *  m  i  1 1 
psi  -  nn*t 

c  *  *  **************** 
c  wi i te ( *, DO )  nn,ps i 

c  DC  formal  {  lx,  ’  nii  =  ',elD.7,'  psi  -  *,elD.7) 

******************* 

c  d  o  8  0  0  i  =  1 , 1 1 

r  write(*,90C)i,ro(i,l),ro(i,2),ro(i/3) 

c90G  fornial(lx,i4,3(el5.7)) 
c  fc  0  u  continue 

c  do  8  D 0  i - 1 ,  n 

c  write(*,900)i,vo(i,l),vo(i,2),vo(i,3) 

c  8  D  0  continue 

********************* 


E  V ALU  AT I  NO 

T 

HE  FOUR  FHI  MATRICIE3. 

phirr 

(1,1) 

= 

4.-3. 0  *cos (psi) 

phirr 

(1,2) 

- 

C  .  0 

piii  r  r 

•  3,3} 

- 

i  .  0 

phirr 

(2,1) 

= 

(  i  .  C*  ( s  in  (  ps  i  )  -ps  i  )  ) 

phirr 

(2,2) 

1  .  0 

phirr 

(2,3) 

- 

c .  c 

phirr 

(3,1) 

= 

(j  .  0 

pli  i  r  r 

(3,2) 

- 

0  .  0 

phirr 

(3,3) 

- 

c  o  s  (  p  s  i  ) 

ph  i  r  v 

(1,1) 

- 

s  i  ri  (  ps  i  )  /nn 

p.h  i  r  v 

(1,2) 

- 

(  2 . 0*  ( 1  -  0  -  cos  (  ps  i  )  )  /nii ) 

l  h  i  r  v 

(1,3) 

- 

0.0 

p  h  i  r  v 

(2,1) 

- 

( 2 . 0  * ( cos ( ps i ) - 1 . 0 ' /nn ) 

ph  i  r  v 

»  z  ,  2  ) 

- 

4.0/nn*sin(psi)-(3.0/nr.*psi 

ph  i  r  v 

U ,  3  ) 

= 

0 . 0 

ph  i  r  v 

(3,1) 

- 

0  .  U 

ph  i  r  v 

(  3 ,  -  ) 

- 

0.0 

ph  i  r  v 

(3,3) 

- 

t.  i  1 1  (  p  s  i  )  /  n  n 

ph  i  vr 

il/D 

- 

3 . 0*  nn  *s i n ! ps i  ) 

ph  i  vr 

(1,2) 

- 

0.0 

p  h  i  v  r 

(1,3) 

- 

0  .  c 

ph  i  vr 

(2,1) 

= 

(  G  .  0  *nn  *  (  cos  ps  i  )  - 1  .  )  ) 

ph  i  vi 

(2,2) 

- 

0 . 0 

p  h  i  v  r 

(2,3) 

- 

0 . 0 

ph  i  vr 

(3,1) 

= 

t  .  0 

p  h  i  v  r 

(3,2) 

- 

0 . 0 

ph  i  vr 

i  3  ,  3  j 

=■ 

-nn*s  i  n  (  p:  i  , 

pi.  i  VV 

(1,1) 

- 

con ' ps i  ) 

/  -3 


c 

c 


phivv(l, 2) 
ph i vv ( 1 , 3 ) 
phi vv ( 2 , 1  ) 
phi vv ( 2, 2 ) 
ph i vv ( 2 , 3  ) 
phi vv ( 3 , 1 ) 
ph i vv ( 3 , 2  ) 
phivv( 3, 3  ) 


2 . *s i n ( ps  i  ) 
0  .  0 


( ~  2 . *s in ( ps  i  )  ) 

-  3 . +  4 . 0 *cos ( ps  i  ) 
0 . 0 
0 . 0 


0 . 0 


cos ( ps  i  ) 


MOVING 

R ( i )  and  V ( i )  FORWARD  IN 

TI 

MF 

do  100 

i  =1 ,  n 

s  unt  1  =  C  . 

0 

s  u  rn  2  =  C  . 

0 

s  u  m  3  =  C  . 

0 

s  u  nr.  4  -  0  . 

0 

s  u  ri  5  -  0  . 

0 

sumO -  0 . 

0 

do  ISO 

3-1/  3 

surnl  - 

suml  +  phir x ( 1 , j  )  *ro (  i  ,  j ) 

+ 

ph i I V ( 1 , 

sum2  - 

sum2  +  pli  i  r  r  (  2 ,  j  )  *r  o  (  i ,  3  ) 

+ 

p  h  i  r  v  1  2 , 

surd  = 

surti3  +  ph  i  r  r  (  3,  j  )  *r.  o  (  i  ,  j  ) 

ph i r  v ( 3, 

sum4  * 

sum*?  +  phivi  (1,  j  )*ro(  i ,  j  ) 

4 

ph i vv ( 1 , 

s  u  m  S  - 

surr.5  i  ph  i  vr  (  2 ,  j  )  *r  o  (  i ,  j  ) 

+ 

phi vv ( 2, 

oUiiiC  - 

suiT.C  i  phivr  (  2,  j  )  *r  o  (  i ,  j  ) 

+ 

ph i vv ( 3 , 

ISO  corititiue 

SUIIi-i 

t  (  i ,  2  ) 

=  s  u  rn  2 

r  (  i  ,  3  ) 

=  sum3 

v(  i,  1  ) 

-  s  u  m  4 

v  (  i ,  2  ) 

=  S  U  fTi  S 

v  (  i ,  3  i 

=  suniC 

10  0  continue 
return 
end 
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j  }  *  v  o  (  i  ,  j  ) 
i ) *vo ( i ,  j ) 
j  )  *vo  (  i  ,  -1 ' 

i )  *vo l  i ,  j  , 

j )  *vo  f  i  , j  ) 
i  )  *  v  o  (  i ,  j  ) 
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phased  arrays  were  rectangles  and  disks,  and  the  other 
arrays  Included  cones  and  spheres.  The  number  of  emitters 
in  the  configurations  ranged  from  five  to  two  hundred. 

The  beam  quality  parameter  was  the  half-power 
beamwidth  as  determined  in  discrete  directions  throughout 
the  hemisphere  below  the  array.  An  array  was  considered 
feasible  if  the  beamwidths  degraded  only  slightly  during  an 
orbit. 

The  most  useful  configuration  was  found  to  be  a  sphere 
of  randomly  spaced  emitters  because  it  provided  narrow 
beamwidths  in  all  directions  below  the  array.  Because  of 
the  motion  of  the  array,  the  beamwidths  changed  in  a 
predictable,  periodic  manner.  The  half-power  beamwidths 
actually  improved  in  certain  directions  and  only  slightly 
degraded  in  others.  The  sphere  also  required  the  least 
number  of  emitters  and  could  be  used  at  any  altitude. 
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